Cosmological Scaling Solutions with Tachyon: Modified Gravity model 
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Modifying the Einstein's gravity at large distance scales is one of the interesting proposals to 
explain the late time acceleration of the universe. In this paper, we analyse scaling solutions in 
modified gravity models where the universe is sourced by a background matter fluid together with 
a tachyon type scalar field. We describe a general prescription to calculate the scaling potential 
in such models. Later on we consider specific examples of modifications and apply our method to 
calculate the scaling potential and the scale factor. Our method can be applied to any modified 
gravity model, in presence of tachyon field. 



I. INTRODUCTION 

Accelerated expansion of the universe has always been 
one of the most challenging subjects of investigation in 
cosmology. For early universe, it plays an important role 
for solving major problems in standard cosmology like 
horizon and flatness problem. It also provides an ef- 
fective mechanism for producing a nearly scale-invariant 
spectra for quantum fluctuations of the inflaton which 
can work as a seed for the structure formations in our 
universe. This paradigm which is dubbed as "Inflation" 
plays a major role in modern era of cosmology [1,]. On 
the other hand, there are not many daunting questions in 
cosmology today than what is the driving force behind 
the late time acceleration of the universe [2]. First di- 
rectly observed by Supernova type la measurements by 
two independent groups in 1998 Q , today it has also been 
confirmed by observations from new SupernovaType la 
measurements|3. Cosmic Microwave Background Radi- 
ations (CMBR) 5|, and galaxy clustering by large scale 
redshift surveys [6 1. 

To discover the nature of driving force behind this late 
time acceleration of the universe, also termed in literature 
as "dark energy" , is one of the most important tasks in 
cosmology today. The first and the most simple choice 
for dark energy is the vacuum energy or cosmological 
constant A. However this possibility of A being dominant 
component of the total energy density of the universe, 
has problem that the energy scale involved is lower than 
normal energy scale of most particle physics models by a 
factor of ~ 10-123. 

An optimistic choice for modelling the missing energy 
of the universe is in the form of a scalar field with a canon- 
ical kinetic energy, slowly rolling down a considerably flat 
potential so that its slowly varying energy density mimics 
an effective cosmological constant. This form of missing 
energy density is popularly called " Quintessence" [7,] . It 
is similar to the inflaton field with the difference that it 
evolves in a much lower energy scale. Amongst them, 
tracker quintessence models have the advantage of allow- 



ing the late time accelerating epoch to be reached from 
a wide range of initial conditions. In order to have a vi- 
able dark energy model, it is important that this scalar 
field remains subdominant during most of the thermal 
history of the universe. In this sense, it is necessary that 
this scalar field mimicks the background energy density 
in the early time. This property is called the "scaling 
or attractor property" . This property is very important 
in cosmology as they allow us to study the asymptotic 
behavior of a particular cosmology and check whether 
that particular behavior is stable or not. There have 
been numerous investigations involving scaling solutions 
in cosmology in standard Einstein's gravity as well as in 
scalar tensor gravity 

Recently, there are suggestions that the present ac- 
celeration of the universe is not due to any new un- 
known component in the cosmic soup, but due to the 
modifications of the gravitational physics at large dis- 
tance scales. In this regard, possible modifications in the 
Friedmann's equation have been proposed. Most of these 
modifications have been inspired by the higher dimen- 
sional brane- world models. A few suggestions are partic- 
ularly interesting in this regard. The Dvali-Gabadadze- 
Poratti (DGP) brane induced gravity modelllo|. the Car- 
dassian model proposed by Freese and Lewis [ll|, a model 
by Dvali and Turner [III, by Shtanov and Sahni[ll| and 
also the modified Chaplygin gas model by Barriero and 
Sen[l3| are some of them. Cosmological solutions with a 
canonical scalar field have been studied in the modified 
gravity models 0]. 

Recently, an effective scalar field theory governed by 
a Lagrangian density with a non-canonical kinetic en- 
ergy term C = -V{4i)F{X), where X = -(l/2)(9''0a^0, 
has attracted considerable attention in cosmology. Such 
model can lead to a late-time acceleration and is called 
"k-essence" [l^ . A scalar field with a non-canonical ki- 
netic term has also been investigated for an early universe 
inflationary scenario and is termed as "k-inflation" [l^ . 
One example of such Lagrangian density is Ctach — 
'V{(j))y/1 — As discussed by Padmanabhan 
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and Roy Choudhury[l7|, this is generalization of the 
Lagrangian of a relativistic particle. The Hamiltonian 
structure of tachyonic matter given by the above La- 
grangian density is very similar to that of a special rel- 

7^ where m 



ativistic particle governed by L 



2 



and q are the mass and generalized coordinate of the par- 
ticle respectively. This tachyon field can naturally arise 
in open string theory [l^ and can provide a rich gamut 
of possibilities in cosmological context [T^. Padmanab- 
han and choudhury ^17j also explored the possiblity of 
this kind field acting both as a clustered dark matter 
and smooth dark energy with a scale dependent equa- 
tion of state (see 20] for a different approach for such 
possibility) . 

In this paper we have studied systematically the scal- 
ing solutions involing a tachyon field described by the 
Lagrangian density,mentioned above, in modified version 
of Friedman equation. In section 2, we describe the equa- 
tions of motion in the modified version of gravity in gen- 
eral and define the variables that allow us to study the 
scaling solutions. The conditions for scaling behavior 
are introduced in Section 3 and we also describe how 
to calculate the potentials for the tachyon field which de- 
scribes the scaling behavior. In section 4, we apply our 
method to some particular examples of the modified grav- 
ity e.g the Cardassian model, the Shtanov-Sahni model, 
the Randall-Sundrum model and the DGP model. We 
conclude in section 5. 



constant. Then the total energy density of the universe 
p, appearing the equation (1) is given hy p — + p^. 
One can calculate the energy density and pressure for 
the tachyon field (p from the Lagrangian density given in 
equation (2) as 



(3) 



We also assume that as in standard gravity, the energy 
momentum tensors of matter and dark energy are covari- 
antly conserved separately which implies that 



—3Hjpj 



) + y(l-'^') = 0. (4) 



Equations (1) and (4) close the system of equations. We 
now define three new variables: 



II. EQUATION OF MOTION IN MODIFIED 
GRAVITY 



We consider a spatially fiat Friedmann-Robertson- 
Walker (FRW) model. In this case the expansion rate 
of the universe H (also called the "Hubble parameter") 
is defined as 



= [8T:G/i)pL\p), 



(1) 



where H = - is the Hubble parameter, a is the scale 
factor, p is the total energy density of the universe. A 
"dot" denotes derivative with respect to cosmic time and 
G is the gravitational constant. Whatever modification 
to standard Einstein's gravity we assume, is parametrized 
by the correction term L(p) and this is assumed to be pos- 
itive definite without any loss of generality. For L{p) ~ 1, 
one recovers the standard Einstein's gravity. 

Our aim is to investigate models where the universe is 
sourced by the matter component together with a dark 
energy part. In this paper, we assume that this dark 
energy field has a non-canonical kinetic energyof Dirac- 
Born-Infeld form with a Lagrangian density 



ach 



(2) 



where V{(j)) is the potential for the field 0. To start with, 
we assume the matter part is described by a barotropic 
fiuid with a equation of state Pj = (7 — 1)P7, 7 being a 



X = 








Y = 








TV = 


Log[a 



(5) 



When we expressed in terms of these new variables, the 
system of equations (1) and (4) becomes 



X' = -{I- X^)[iX - y/WX] 

V3A2xr(r - 3/2) 



- 37 



(6) 



A 



3A 



7 - 



dLog\L{p)\ 
dp ' 



where "prime" means differentiation with respect to N . 
The parameter A and F are defined as 



A 
F 



{dV/dcj)) 



VS^L{p)Vy^ 
Vd^V/d(P^ 
(dV/dc/))^ ■ 



(7) 



From the definition of total energy density, one can get 
the constraint equation 



Y2 



p 



1. 



(8) 
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The equation of state for the tachyon field and its den- 
sity parameter are given by 



70 = X' 



VLa, — 



y2 



VI 



(9) 



Then the allowed phase space for X and Y is given by 
< X"^ -|- < 1 from the requirement < fl^ < 1. 
From the expression of the equation of state 7^, one see 
that 70 > (F^ = (70 - 

For the case A ~ constant^ in standard gravity, the 
potential turns out to be V{<t)) oc 0^^. But due to the 
presence of the modified term L{p) ^ 1, it is not so in our 
case. But it is interesting to see that, for A = constant^ 
the equations in (6) have the identical form to that of the 
plane autonomous system of standard relativistic cosmol- 
ogy involving a tachyon field in terms of the parameters 



Xs = 
Ys = 



V87rGV^((/)) 
{dV/d(j)) 



(10) 



This immediately ensures that our system of equations 
(6) with constant A admits the same set of critical points 
to that of the standard cosmology when those solution 
are expressed in terms of {X, F, A}. Also the stability of 
these critical points can be directly obtained from sta- 
bility analysis of the standard cosmology scenario. Now 
there are altogether five fixed points in our system equa- 
tions (6) when {X,Y,X] = {Xc,Yc, Xc} arc constants. 
Three of these {X^ = 0,Yc = 0}, {X^ = 1,^^ = 0} and 
{Xc = —l,Yc = 0} are unstable as studied by Copeland 
et al[2l[. The fourth fixed point is given by 



Xc 
Yr 



7, 

3"7/Ae 



(11) 



^ yA4T36-A? 



This fixed point is a stable node for 7 > a 
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III. SCALING SOLUTIONS 

As we mentioned in the earlier section, that for a stan- 
dard cosmology, A = const, corresponds to a potential of 
inverse power law form. But with the modification term 
L{4>) in the expression of A in equation (6), the form of 
the potential will change. In this section, we shall cal- 
culate the form of the scaling potential in the modified 
gravity scenario. For this purpose, we essentially follow 
the method earlier prescribed by Copeland et aZ.[22] for 
a standard canonical scalar field. 

One can show that for both the critical points (11) and 
(12), X' = r' = a' = 0, if the condition 



T=^-+p^^Log[L{p)] 



(13) 



is satisfied. Using the relation p = V/Yc, the above 
equation can be written as 

This is very similar to the corresponding equation (eqn 
(19) in ref [231) ^ canonical scalar field except the fac- 
tor (3/2) in the second term in left hand side of the above 
equation. The above equation can be integrated to find 
p{(j)) which together with the expression V = Y^p will 
give the corresponding potential for the scaling solution 
for given choice of modification L{p). By integrating the 
above equation and using equation (7) with constant A, 
we get 



dp 



-nXrYr' 



(15) 



It is a stable one for < 7 < a = ^ yJXc + 36 — A^ 

For this one can easily show that 7^ = 7, hence it rep- 
resents a scaling solution. But in this case as Vl^ — 
37/ (A^Vl ~ 7); 7 h^-s to be less than 1. It restricts this 
case to be a viable model for dark energy as the back- 
ground fluid can never be matter-like. The fifth fixed 
point is given by 



Xr 



Yc 



Xc / y/AfTse-Ag 

V3 I 6 



1/2 



VAf+36- A^ 



1/2 



(12) 



where — 8ttG. While deriving the above expression, 
we have put one integration constant to zero without no 
loss of generality. This is equivalent to giving a linear 
shift to the value of the scalar field. 

It is also interesting to calculate the behavior of the 
scale factor for a given scaling solution assuming the 
scalar field is monotonically varying function of time 
((/) ^ 1). In general one can express the scalar field equa- 
tion (4) in the following manner: 



= (-3iJ0=^ 



-3i702 



(16) 
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By using the definition of the Hubble parameter H = ^, 
one can write 



37^2 ^ 1 da dp^ 1 



(17) 



Substituting the above equation in the Hubble equation 
(1), one can write 



3.0 
2.5 - 
2.0 - 
£ 1.5 
1.0 



da dp 2 2r2^i,\ 

— — = —K ap L (0), 



(18) 



where p is the total energy density of the universe and 
can be written as p = . Defining a new variable 6(0) 



h(4>) — exp 



dp 



p'L^p) 



(19) 
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FIG. 1: The evolution of the universe for Randall-Sundrum 
II brane world models. We have taken = 10"^" for this 

figure. Ac — 1, 1.5, 2 from top to bottom. 



one can write equation (18) as 
da db 



-K ab. 



(20) 



which can be subsequently integrated to give the scale 
factor a: 



a — Exp 



.2/db.i 



(21) 



Using the equation (16), one can also calculate the time 
dependence of the scale factor through the following 
equation: 



t^-^/3K / d0p3/2(0)L(0)(Jl) 



dp. I 



(22) 



IV. DIFFERENT CLASSES OF MODIFIED 
GRAVITY 

A. Randall-Sundrum model 

For Type-II Randall Sundrum brane-world model 
where 3 brane with positive tension is embedded in five- 
dimensional Anti de-Sitter spacetime, the modification 
to the standard gravity is given by 



(23) 



a is the tension of the 3-brane. One can now use equa- 
tions (5) and (15) to determine the scaling potential 
which is given by: 



aX^^Y^^KU^ - 2 



(24) 



One can also obtain the time dependence of the scale 
factor a{t) using equations (21) and (22): 



a{t) 



6 



1 



(25) 



As explained in section H, there are two stable solu- 
tions. One of them can not produce a viable late time 
dark energu model, as the background fluid can never 
behave like matter. For the other case, one can have a 
late time scalar field dominated case {ft^ = 1) with ac- 
celerated expansion. We have plotted the time evolution 
for this case in figure 1. We have taken = 10"^°. We 

have shown the evolution for different values of A. We 
should mention that evolution of the universe is same as 
obtained by Copeland et al. in [22i] for a standard canon- 
ical scalar field. But the due to the noncaninical nature 
of the kinetic term, the required potential is different. In 
this case the acceleration occurs for 7^ < 2/3 in the late 
time when the p term in the action dominates whereas 
for 70 < 1/3, acceleration occurs in early time when the 
p^ term dominates. This translates to the constraints on 
Ac as Ac < 1.86 and Ac < 1.01 respectively. 



B. Shtanov-Sahni Braneworld Cosmology 

There is another interesting class of models proposed 
by Shtanov and Sahni where a 3 brane with negative 



5 



tension a is embedded in a five dimensional conformally ! 
flat space, where the fifth dimension has time-hke signa- 
ture. In this case the modification to standard gravity is ^ 
describe as: 



One can calculate the potential which is given by 



(26) 3 - 



(27) 



One can also calculate the corresponding scalar factor 
a{t) 



(28) 
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FIG. 2: The evolution of the universe for Cardassian model. 
We have taken -A: — 10^^" for this fieure. Ac = 1, 1.5, 2 from 
top to bottom. 



C. Cardassian model 

All the above modified gravity models modifies the 
gravity at small distance scales i.e it modifies gravity at 
high energy scales (early times). But recently Snia, CMB 
and LSS measurements have confirmed the late time ac- 
celerated expansion of the universe. A number of phe- 
nomenological models have been proposed to explain this 
late time accelerated expansion through modifications of 
standard gravity at large distance scales. In one such 
model, know as "Cardassian model" originally proposed 
by Freese and Lewis the modification is given by 



L{p) = + 



(29) 



The present acceleration of the universe can be obtained 
in such model with n < —1/3 when the universe contains 
only the matter. Although the presence of any scalar 
field was not assumed in this original model, it is inter- 
esting to study the scaling behavior of the background 
cosmology in the presence of both matter and a scalar 
field, a noncanonical one for our present purpose. To get 
the result analytically in closed form we have assumed 
n — —1/2 in the following calculations. We should men- 
tion that taking different values of n is always possible, 
only problem is that solution may not be in a closed form. 
We have used this particular value of n to present the so- 
lution in a closed form. With this choice of n, one can 
directly obtain the form of the potential F(0) and the 
corresponding scale factor a{t) as following: 



F(0) 



(30) 



a{t) 



A 



4\/3 



(31) 



We have shown the evolution of universe for this in 
figure 2. We have assumed A/m4 = 10^^". Acceleration 
of the universe takes place in early time for 7 < 2/3, 
which constrains the Ac parameter as Ac < 1.86. For late 
time, universe always accelerates. 



D. Dvali-Gabadadze-Porrati Brane-World Gravity 
Model 



The Dvali-Gabadadze-Porrati (DGP) model is a brane- 
world model (a 3-brane embedded in 5D anti-deSitter 
bulk) with the incorporation of the scalar curvature term 
for the 3-brane in the total action. The bulk is empty and 
all kind of matter fields are restricted to stay only on the 
3-brane. The modified FRW equation for DGP model is 
given by 



ro 



(32) 



where Tq 



Tj-T- with 7714 and ms being the 4-d and 5-d 

Planck mass respectively. The parameter ro actually de- 
termines the scale at which the switchover from standard 
gravity to modified one takes place. The (-I-) and (-) sign 
corresponds to different kinds of embedding the 3-brane 
in the 5d bulk. The modified equation now becomes 



H 



1 

2^ 



Tl + Vl + aiP 



(33) 
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where ai = A direct comparison of the above 

equation with equation (1), results the correction term 



in the Friedmann equation as 



L 



1 



aip 



(34) 



Putting this in the equation (15), one can integrate to 
get 



=Fl + VI + "iP 



aiSinh ^ 



1 



= KX,Y,(j) (35) 



This equation together with equation (5) can be used to 
find the scahng potential V{(t)). One can also calculated 
the corresponding time dependence. This is given by 



3X| 
2ro 



t = Cotanh ^■\/l + aip + 



1 



(36) 



V. CONCLUSION 

The issue of late time acceleration of the universe has 
been one of the most serious challenges in cosmology 
nowadays. While including an extra dark energy compo- 
nent with repulsive gravity in the energy budget of the 
universe, is one of the most studied approaches for ex- 
plaining, modifying Einstein's gravity at large distance 
scales, has also been taken seriously in recent times. 
While any modifications to Einstein's gravity has its own 
problems, this idea has been rigorously pursued by var- 
ious researchers in recent times. One of the main moti- 
vations of such idea is that, modifications of Einstein's 
gravity can arise at different energy scales very naturally 
through compactifying higher dimensional theories. 



On the other hand, the concept of scaling solutions 
in cosmology has been taken seriously in recent times, 
as such solutions are necessary for solving cosmic coin- 
cidence problem in dark energy models. In this paper, 
we have studied in systematic way the scaling solutions 
in modified gravity models when the universe contains a 
tachyon type scalar field in addition to the sandard mat- 
ter field. This is an extension of the earlier work done by 
Copeland et al. [l^l where the scalar field with a canonical 
scalar field was considered. 

We first describe the general equations and method of 
calculating scaling potential. Later on, we take specific 
modified gravity models and apply our method. We take 
four choices of modifications, e.g the Randall-Sundrum 
II model, Shtanov-Sahni Model, Cardassian model and 
DGP model. Our method can also be used to con- 
sider any modified gravity to calculate the scaling po- 
tential, while considering the tachyon type scalar field. 
We should mention that Tsujikawa and Sami[23j have 
earlier considered scaling solutions in modified gravity 
with a tachyon field. But for the modification, they 
have considered the special case when cx p" which 
does not include modifications like DGP model. Das et 
al.[2^ have also considered tracking solutions in modi- 
fied gravity both with quintessence and K-essence type 
fields. Their approach is different from ours and they 
have not considered fields like tachyon which has some 
specific features. 

In future, one can take the general K-essence action for 
the non-canonical scalar field, and calculate the scaling 
potential. This will be our future goal. 
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